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Quantum dynamics of in-plane Josephson junctions between two d-wave superconducting films is 
described by the anisotropic XY-mode\ where both, quasiparticle and Cooper pair tunneling terms, 
appear to be equally non-local. Applying a combination of the weak and strong coupling analyses 
to this model, we find a compelling evidence of a dissipative phase transition. The corresponding 
critical behavior is studied and contrasted with that found previously in the conventional (s-wave) 
Josephson junctions. 

Quantum phenomena in ultrasmall Josephson junctions between close pairs of superconductors have long been at 
the forefront of research in theoretical and experimental condensed matter physics citesz. However, putting a few 
exceptions aside, the previous theoretical work was almost exclusively focused on the Josephson junctions between 
conventional s-wave superconductors which are fully gapped. 

The growing interest in unconventional (e.g., d-wave) Josephson junctions has been bolstered by the studies of 
cuprates and other gapless superconductorsS*^ as well as the proposed applications of mesoscopic d-wave devices in 
quantum computing^. 

A unifying theoretical description of tunnel junctions between both, normal metals and superconductors, can be 
constracted in terms of the imaginary-time effective action which governs the dynamics of the phase diference across 
the junction (/)(r) — <f)R{T) — (f)L{T) conjugate to the junction's charge g(r)i. Under the assumption of weak tunnehng, 
this action can be written in the general form 

AE'cJo ^Ih^ ~ X io f"^^^ ~ ''^^ 2 ^^^^ ~ 2 I'^^i'^^a (1) 

where T is the temperature and the first (local) term is the charging energy (measured in units of Ec = e^/2C) of a 
junction with capacitance C, while the (potentially non-local) a- and /3-terms represent the processes of quasiparticle 
and Cooper pair tunneling across the junction, respectively. The corresponding integral kernels 

a{T) = - ^ \tk,k' |'Gl(t, k)GR{~T, k') 

k,k' 

P{t) ^Y.\^k,k'?FL{T,k)FR{~T,k') (2) 

k,k' 

are determined by the Fourier transforms of the normal and anomalous quasiparticle Green functions on the left /right 

bank of the junction (GL,_R(iw„, /c) = (iw„ -t- ^l,_r)/(w2 + ^2 + A|_^) (md FL,R{iuJn,k) ^ Al,r/ {cjI + + Al j^^), 

respectively) and the tunneling matrix element tk,k' ■ 

Under the standard assumption of a generic, momentum-non-conserving, tunneling [tk,k' ~ const] between two fully 

gapped (s-wave) superconductors, the pair tunneling kernel decays exponentially (/3s(t) oc exp(— 2Asr) for r ^ 1/As). 

As a consequence, the constituents of a Cooper pair tunnel almost simultaneously, and the /3-term reduces to the 

\ It \ It 

single time integral —Ej cos 0(r)dT proportional to the local Josephson energy Ej = j3{T)dT. 

Although the same argument may seem to be equally applicable to the quasiparticle tunneling (a-) term, a proper 
description of the experimental data on realistic junctions requires one to retain a non-local qusiparticle tunneling 
term whose kernel behaves as ols{t) oc I/t^ in the entire range l/A^ < r < 1/T. In this way, one can account for a 
finite subgap resistance due to inelastic (phase-breaking) scattering!. 

In contrast, a derivation of the effective action (1) of a junction between two gapless d-wave (as well as any I ^ 0- 
wave) superconductors along the same lines manifests its non-universality and a strong dependence on the details 
of the tunneling matrix elements. For a momentum- independent tunneling the kernel (3{t) given by Eq.(2) vanishes 
identically, alongside the angular averages performed in the course of the momentum integrations Fl^r{t, k). 

By contrast, in the presence of a momentum-conserving [|ife,fc'P ~ S'^{k — k')] node-to-node tunneling across the 
junction between two three-dimensional d-wave superconductors both kernels in Eq.(l) demonstrate the algebraic de- 
cay ad,3D ~ l3d,3D oc I/t"^, thereby resulting in the super-Ohmic quasiparticle dissipative term in Eq.(l)^. Moreover, 
the Ohmic quasiparticle tunneling term might still occur for those relative orientations between the d-wave order 
parameters that allow nodal quasiparticles to tunnel directly into the surface-bound zero energy states^. 

Interestingly enough, inspite of its being of an equally non-local nature, the kernel of the concomitant Cooper 
tunneling term is routinely treated as though it was local, and it has become customary to endow the effective action 
(1) with the ad hoc conventional Josephson energySi^. 
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The latter approximation may indeed prove adequate for the junctions formed by three-dimensional grains of a 
d-wave superconductor, as, e.g., in the case of strongly inhomegeneous high-T^ samples^. However, it would obviously 
fail in the case of tunneling between a pair of genuinely two-dimensional d-wave superconductors (e.g., thin films) 
where both the quasiparticle and Cooper kernels show the same Ohmic decay^ 

ad,2_D=a/T^, l3d,2D^f3/T^ (3) 

In the presence of elastic scattering, this Ohmic behavior of both q(t) and /3{t) (see Erratum in Ref.^) will undergo a 
crossover to the exponential decay at time scales in excess of the inverse bulk impurity scattering rate 7, thus effectively 
restoring the local Josephson energy. However, as suggested by the wealth of transport data in quasi-planar cuprates, 
7 turns out to be quite low compared to the maximum gap A^, and, therefore, the kernels (3) remain essentially 
non-local in the entire regime max[l / A^, l/i?c] < t < min[l/T, I/7]. 

Apart from its anticipated relevance to the problem of tunneling between thin d-wave films, the anisotropic XY- 
model described by the effective action (1,3) is interesting in its own right, for it appears to reside outside the realm 
discussed in the literature up to date and its quantum dynamics has not yet been studied. To fill in the gap, in what 
follows we investigate this new model by applying a number of complementary techniques that cover both weak and 
strong dissipative coupling regimes. 

In the weak coupling regime (a, /3 < 1), an adequate approach to the model (1,3) is provided by a direct perturbative 
expansion for the grand partition function in the charge representation 

rQ+ne 00 N \jT fXjT 



00 M „i/T j-l/T 1 j-l/T 



(r)dr) (4) 



where the lower limit in all the integrals is set at Tc ^ 1/^'c, and the instantaneous value of the total charge of the 
junction 

JV M 

q{r) = Q + ne + e^(0(r - r+) - ^?(r - rr)) + g,[(?(r - r+) + e{r - r")] (5) 

includes a continuously varying contribution Q — CV induced by an applied external bias V . 

The sum in (4) accounts for all the trajectories in the charge space {Q + ne ^ Q + {n± l)e . . . ^ Q + ne) that 
consist of N pairs of quasiparticle {Q ^ Q ± e followed by Q' Q' ^ e) and M pairs of Cooper pair (Q ^ Q ± e 
followed by Q' ^ Q' ± e) tunneling events. 

The periodic dependence of the partition function (4) (hence, any physical observable) upon the external charge Q 
with the period e allows one to restrict its values to the "Brillouin zone" (BZ) — e/2 < Q < e/2 in the charge space. 
Unlike in the case of the local Josephson energy, there is no room for a periodicity with the minimum charge 2e even 
for a ^ 0. 

In the absence of tunneling, the ground and first excited states appear to be degenerate (Q^/2C — (Q±e)^/2C) at 
the BZ boundaries {Q — ±e/2). In the vicinity of the degeneracy points the renormalized gap A{Q) — Ei{Q) — Eq{Q) 
between the ground and first excited states remains small compared to i?c, and, therefore, one can neglect any 
transitions between these two energy levels and the rest of the spectrum separated by an energy gap of order E^- As 
a result, the sum (4) reduces to that over the trajectories comprised of a sequence of AQ = ±e "blips" between the 
two lowest statesZi. 

A straightforward analysis shows that this effective two-state system can be amenable to the renormalization group 
(RG) analysis, akin to that developed in the context of the Kondo and other quantum spin-1/2 problems. In particular, 
the standard procedure of changing the cut-off in the time integrations in (4) from Tc to > Tc and integrating over 
all the pairs of opposite blips with separations Tc < — t~\ < results in the renormalized partition function 
whose relevant part retains the original form (4), albeit with renormalized parameter values. 

To first order in the dissipative couplings, the wave function renormalization Z (defined through the propagators 
of the two lowest states, Go,i — Ze~^°-^'^) and the vertex renormalization factor Zy receive contributions of order 
a In t^/tc and I31tit^/tc, respectively (notably, the closed Dyson-type equations for the propagators Go,i(t) and the 
vertex function F contain no analog of the polarization function—). 

Combining these results, one arrives at the RG equations for the effective couplings Ur — aZ^Z^, (3r — (3Z^Z^, 
and Ar = A^Tc 

— = ~2ar{ar - Pr) 

dmTc 



3 



df3r 

— 2pr(ar — Pr 



dluTc 

d In A, 



d\n: 



(1 - 2ar) (6) 



which are vaUd for as long as the renormaHzed dissipative couphngs and Pr remain small. 

By solving the RG equations (6) and evaluating their solutions at the lowered energy cutoff l/r^ = one obtains 
the renormalized couplings 

/3 l + 2?;ln£;c/A^ ^ ' 

where rj = a — p. In addition, the RG procedure yields a self-consistent equation for the renormalized gap 

" [l + 277ln£:e/A^]"/'' 

the bare value of which is Ao(Q) = e[°\q) - E^°\q) = E^l - 2Q/e). 

As follows from (7), for > the RG trajectory flows towards weak coupling (see Fig.l) where the invariant 
charge A^ increases, although the actual gap A^ given by (8) continues to decrease. In this regime, the quantum 
current fluctuations associated with the phnomenon of Coulomb blockade completely destroy the classical Josephson 
effect, and the junction remains in the insulating state. This behavior appears to be similar to that emerging in 
the previously studied limit /? = where the system possesses the exact XK-symmetry and remains insulating for 
arbitrary values of the external charge Q-^. 

Near the line rj = the RG flow slows down, and right on that line the effective coupling — Pr undergoes 
(almost) no rcnormalization, while the effective splitting demonstrates a power-law dependence on its bare value 

A, = Ao(Ao/£;c)'"/i-'", (9) 

which behavior is indicative of a possible dissipative phase transition at a = P = 1/2. 

Exactly at = the system acquires the Ising symmetry, which suggests that the transition in question is likely to 
be of the Kosterlitz-Thouless typc^. Moreover, since the effective coupling a — P undergoes no rcnormalization (at 
least, to first order), this critical behavior appears to be reminiscent of that of a junction with the local Josephson 
energy and Ohmic (as opposed to the quasiparticle tunneling-induced) dissipation. 

In the regime a = P < 1/2, the quantum phase fluctuations are strong and the energy bands i?o,i('3) remain non- 
degenerate, thereby forcing the junction to operate in the insulating (Coulomb blockade-dominated) regime. However, 
for a = P > 1/2 the quantum fluctuations get quenched, and the energy bands become progressively more and more 
degenerate in a finite portion of the BZ which expands from the boundaries Q — ±e/2 inward as the parameter a — P 
increases. This behavior indicates a suppression of the Coulomb blockade and a possible restoration of the classical 
Josephson effect where the voltage drop across the junction vanishes and the current is determined by a (nonzero) 
average value of the phase difference < 0(t) >= const of the superconducting order parameters on the opposite banks 
of the junction. 

In contrast, for rj < one finds a runaway RG flow towards strong coupling (see Fig.l) where both and Pr 
become of order unity and A^ starts to decrease, regardless of the bare values of these dissipative couplings. Such a 
behavior suggests that the junction is likely to remain in the (super) conducting regime for all a < p. 

Using the formula (8) for the gap, one can readily compute the ground state energy Eo{Q) = Ec — ^Ar{Q) which, 
in turn, yields the junction's effective capacitance 

Cr = i^^^r'lQ^o - C[l + Arj + Oiv^r/^ (10) 

Even more pronounced, appears to be the effect of tunneling on the average excess charge of the junction close to the 
BZ boundary {5Q = Q - e/2 0") 

/ N ^ ^dEciiQ) e, siqnSQ . , , 

< '^(^^ ^ - - 2^' + [l + 2,hr(e/^g)]^A> ) ^''^ 



Consistent with the anticipated dissipation-induced suppression of the Coulomb blockade, Eq.(ll) demonstrates the 
rounding of the classic Coulomb staircase-like dependence of < q{T) > on Q with increasing a. 
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Also, by continuing Eq.(8) analytically to the regime Q > e/2 where the logarithmic function of 5Q acquires an 
imaginary part, one can estimate the decay rate of the first excited state with energy Ei{e/2 — 6Q) — -Eo(e/2 + 6Q) 

27ra(SQ/e) , , 

r(Q) . ImE^iQ) . i^c ^^^^^^^^^^^,^^)/,,,,,,, (12) 

The excited state remains well defined for as long as T{Q) <C Ar{Q), which condition can only be satisfied in the 
weak coupling regime (77 ln(e/i5(5) < !)• 

Notably, the effects of the external charge screening and the reduction of the tunneling rate manifested by Eqs.(ll) 
and (12), respectively, are both enhanced as compared to the case of a normal junction ("single electron box") where 
/3 = 0l. 

These two effects are also exhibited by the current- voltage (/ — V) characteristics. In a voltage-biased junction, the 
induced DC current 

HV) > rmo-vc « - EM (13) 

vanishes below the threshold Vc — Ec/e. However, at finite temperatures the hard Coulomb gap gets partially filled 
with thermally excited quasiparticle excitations, thus giving rise to a temperature-dependent conductance. By the 
same token, in a current-biased junction the zero-temperature I — V characteristics are expected to become non-linear 
for / < eEc. 

Next, we consider the regime of strong dissipative couplings (77 > 1) where the first insight into the problem can 
be obtained by virtue of an adaptation of the variational technique. In this method, the correlation function of the 
small (" spin wave-like" ) phase fluctuations is sought out in the form 

The variational parameters g and D are then used to minimize the free energy F = — TlnZo + T < S — S^Sq > where 
5'^So{(t)u) is the quadratic action which corresponds to Eq.(14) and determines the averages < e"*^ >= exp(— i < 0^ >). 
In this way, one obtains a self-consistent equation for the effective dissipative coupling 

,y = g(9+l/9-l) (2/3)1/1-9 (I5) 

It can be shown that if the following condition is fulfilled 

a - /3 > 60(1/2-/3) In ^ (16) 

[here e is the Euler's number], there exists a finite gap due to the Xy-anisotropy in the two-dimensional vector space 
spanned by the unit vector n = (cos^/2,sin(/)/2) 

D ^ Ec{2(3y/s-^g^/^-s (17) 

where g « ri^^^^'^{2(3y^'^ is the solution of Eq.(15) for rj > 0. For comparison, within almost the entire domain rj < 
one finds g ~ rj. 

The gap (17) attains its bare value 2/3i?c at g — > oo and decreases upon lowering the energy cutoff. Considering 
that the bandwidth of the phase fluctuations, too, undergoes a downward renormalization from its bare value Ec (see 
Eq.(28) below), one finds that the spectrum of the quadratic operator S'^So/ScjP remains stable for both positive and 
negative yy. 

It should also be noted that, when continued to real frequencies, the spectrum of the small phase fluctuations 
appears to be strongly overdamped. Therefore, except for the immediate vicinity of the separatrix rj = (hence, in 
almost the entire a — /3 plane), it bears no resemblance to the simple plasmon pole (cf. Ref.®). 

In the presence of a nonzero XF-anisotropy D ^ 0, the expectation value < (j>^{T) > becomes finite and the phase 
(j) gets localized. As a result, the real part of the AC conductance 

G{u;) = I{u;)/V{lj) = Re ] (18) 

develops a coherent peak at zero frequency, G{uj 0) = DS{uj), thus indicating the onset of the classical Josephson 
effect. 
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In the case that Eq.(15) features no solution, the gap D vanishes and the expectation value of the phase fluctuations 
< ^^(r) > diverges, so that the phase remains delocalized. At the first sight, it may seem that the real part of the 
conductance approaches a finite value G'(O) = 5 in the DC limit. However, we expect that the continuing downward 
renormalization of the effective dissipative parameter 77 (see Eq.(19) below) will eventually bring the system into the 
insulating regime. 

In order to justify the above claim, we point out that is in the absence of a non-trivial mean field solution the 
variational method must be abandoned in favor of a straightforward perturbative approach. The latter yields a 
continuous renormalization of the dissipative couplings due to the non-Gaussian (quartic and higher order) terms in 
the action (1). 

To the leading order, such a renormalization is described by the RG equations 

dar o-r d(3 (3r ^^g^ 



which demonstrate a monotonic decrease of the effective dissipative parameter as a function of the lowered energy 
cutoff 



^^ = r^--l^lnE,T'^ (20) 



It is worth noting that the steady decrease of both ar and Pr as well as the approximate constancy of the ratio ar/ (3r 
manifested by Eqs.(19) are consistent with the behavior found in the weak coupling regime, and so the RG trajectories 
interpolate smoothly between the weak and strong coupling regimes. 

Interestingly enough, the RG equations (19) coincide with those of the spin-boson model 



S'sf,[m(T)] = ^ j dT{l-mz)^+ J dT{^Ecml + m^hx + ruyhy) 



(21) 



where m — {^/T^^^rn^ cos 4>/2, — m\ sin 0/2, m^) is a 3D unit vector field composed of the phase and a conjugate 
"momentum" m^. The first term in (21) is the spin's Berry phase, and a two-component random field hx^yir) represents 
a dissipative bath governed by the correlators < |/i2:(a;)p >= {a + (3)\uj\ and < \hy{uj)\'^ >= r]\uj\. 

The origin of the equivalence between the strong coupling regimes of the actions (1,3) and (21) can be traced back 
to the observation that at large 77 both, the discreteness of charge and the quantum nature of spin (or, for that matter, 
the spin's Berry phase), turn out to be largely irrelevant and the actions (1,3) and (21) are both dominated by the 
(identical) dissipative terms. 

As another side note, we mention that the action (21) can also be encountered in generic (spin anisotropic) models 
of quantum impurities as well as those of noisy qubits (generalized (pseudo)spin-l/2 systems exposed to dissipative 
environments). In light of the earlier observations^ that in such models the only stable fixed points tend to be those 
of the Ising symmetry, it does not seem unreasonable to surmise that the behavior in the entire region defined by 
Eq.(16) can, in fact, be similar to that found for a — (3 > 1/2. 

Altogether, the above results suggest a tentative layout of the phase diagram where the insulating behavior sets 
in throughout the domain defined by the inequalities < f3 < 1/2 and P < a < f] + eln(l/2/3) (see Fig.l). By 
implication, we then conjecture that upon increasing the value of 77 > or changing its sign, the insulator gives way 
to the (super)conducting phase. 

In order to further refine our understanding of the strong coupling regime we employ the semiclassical approach in 
the phase representation*. The Coulomb energy (which, contrary to the weak coupling case, is now subdominant to 
the dissipative terms) induces a splitting between the (otherwise degenerate) phase configurations ("vacua"), and the 
partition function 



Z(Q) '^'^0 / i?0(r)e-^['^(")l+2-^^'3/^ (22) 



is now saturated by the trains of 47r-phase slips (" instantons" ) corresponding to the transitions between different 
even/odd vacua 4>-2m — 2nM and (j)2M+i = (2M -|- l)n. 

A viable candidate to the role of the trajectory connecting the vacua 4>2M{2M+i) and (I)2M+n{2M+i+n) can be 
chosen in the form (here = ±1 and X]fe=i Ik = N) 

n 

A0Ar(r) =4^gfetan-if]fe(r-Tfc) (23) 
fc=i 
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which is inspired by the exact solution of the classical equations of motion for (3 = and Ec = oo. 
Evaluating the action (1,3) on the single instanton ("bounce") trajectory A(/)i(t) we obtain 

5i(f^) = 2A + — + — ^ (24) 

Minimizing Eq.(24) with respect to the instanton's size we find ^min ^ E^fi^l'^ which is small compared to the 
average spacing between the instantons (|t„ — rm| ~ 1/NT) for T — > 0. 

Since the instantons are well separated, one can expand about a single-instanton configuration (0(r) = A0i(r) + 
(5(/>(t)) and obtain a quadratic action of the Gaussian fluctuations 

6^Si = r^(a|(50^n|cj + n/2\ + \uj- n/2\ - + l3[5dp^5^},^SEc/T -\uj + - n/2) 

UJ 

+Sc^^S(l>},+jEjT -luj- n/2\ - n/2)]) (25) 

At /3 = the eigenvalues of the operator with the frequencies |a;| < turn out to be independent of a. However, 
for any finite (3 these "soft" modes acquire gaps and, therefore, can not be separated from the rest of the spectrum 
(cfi^). The fact that this spectrum remains positively defined and, therefore, stable justifies the use of Eq.(24) as a 
sensible approximation for the (unknown) true local minimum of the action (1,3). 

Next, by singling out the zero modes corresponding to the changes in the instanton position {S(j)'{'^ cx d(f)i/dT) and 
size {S(j)\ oc d(t)i/dQ) and integrating over the non-zero normal modes, one can cast the partition function as that of 
a dilute instanton gas 

where J(0) is the measure of integration over the zero modes, and the determinant stemming from the integration 
over the non-zero modes is normalized with respect to that of the operator S'^Sq/S^)^ corresponding to the Gaussian 
fluctuations about the trivial vacuum 0(t) = 0. 

At large positive values of -q the ground state energy can be obtained by taking the zero temperature limit of the 
logarithm of Eq.(26) 

Eo{Q) = -TlnZmr^o ^ E,cos{2^Q/e) j dl7J(f^)(^^^|J|^)'/'e-^^(") (27) 

At Q — Eq.(27) yields a renormalized bandwidth of the phase fluctuations whose asymptotic behavior for 7/^1 
reads 

E,^r - E.r^e-^^"" (28) 

By analogy with the /? = case, this reduction of the bandwidth can be attributed solely to the exponential blow-up of 
the effective capacitance Cr ^ C exp(2Tr^r])/r]. For comparison, along the separatrix a = f3, the effective capacitance 
grows only as Cr ~ Cexp(7r-\/a/8)/a^/^ for large a. 

The advantage of working in the phase representation is that one can better elucidate the nature of the putative 
dissipative phase transition by relating it to a spontaneous breaking of the symmetry between the (j)even and (podd 
vacuaS. 

In the disordered phase, the correlation function < e^^^'^^/^e""^'"'/^ > decays algebraically and the insulator-like 
I ~V characteristics show the presence of a hard (at T = and / — > 0) Coulomb gap. By contrast, in the ordered 
phase the system develops an order parameter < e*"^'"^'/^ >^ whose presence indicates that the phase variable is 
localized in either even or odd vacua. This does not, however, constitute a complete phase localization, and, therefore, 
the apparent classic Josephson effect should be considered a non-equilibrium phenomenon which can only be observed 
at flnite (albeit potentially quite long) observation times, while at still longer times the response would revert to the 
resistive behavior^. 

Another potentially important non-equilibrium effect is excitonic enhancement of the tunneling probability which 
modifies the exponent in the power-law in Eqs.(3) to 2 — e where e is a function of the (non- universal) scattering 
phase shift^°. In this "sub-Ohmic" case, the right-hand-sides of the RG equations (19) pick up additional terms 
ea and e/3. Consequently, one now finds a fixed point at rjc = l/27r^e, upon approaching which the bandwidth 
vanishes as i?* ~ Ec{l — r]/r]cY^~'^^^'^. One would then be led to the conclusion that in the domain defined as 
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T]c < ct — P < eln(l/2/3) the transition from the insulating to the Josephson phase is preempted by that into a new 
conducting (albeit, possibly, dissipative) state. 

We anticipate that the prediction of the dissipative phase transition made in this work might be possible to test 
by using junctions between thin single-crystal cuprate films. For the parameter values found in Ref.® (a — 2(3 and 
Ec ~ Ad) one would expect to observe the apparent superconducting response at a > etc ~ 0.85. Considering that 
the typical junctions manufactured from the cuprate superconductors tend to have relatively high conductances^, this 
condition can be readily met. 

A still more realistic model of the cuprate films would also feature a renormalized kernel q;(t) that accounts for 
possible zero-energy bound states (which is only absent in the case of the node-to-node tunneling) as well as a quadratic 
Ohmic dissipative term representing the effect of a shunting normal resistance (metallic environment surrounding the 
superconducting grains). We anticipate that by including these effects into the model (1,3) the borderline between the 
(super)conducting and insulating phases can be pushed even further towards lower values of the dissipative couplings. 

In contrast to our findings, the naive gradient expansion employed in Ref^ is not capable of revealing any critical 
coupling ac or the corresponding phase transition in question. In fact, a direct application of this technique to the 
model (1,2) would yield an undamped correlation function < \4)u\'^ >— {uj'^/E* + D*)^^ with a strongly temperature- 
dependent gap A ~ (EcTY^'^ throughout the entire a — f3 planeS^, which behavior is starkly different from that found 
in this work. 

Also, as far as a proper modeling of the highly inhomogeneous cuprate samples studied in Refs.^ is concerned, 
the latter are likely to be best described not by the strictly two-dimensional action (1,3) but rather as a network of 
three-dimensional grains whose behavior is expected to be somewhat more conventionali. 

To conclude, in this work we carried out a comprehensive analysis of the quantum dynamics of a node-to-node 
Josephson junction between two planar d-wave superconductors. Our results reveal the presence of a dissipative 
phase transition in the a ~ [3 plane. We also discussed properties of the physical observables that are indicative of 
the transition in question, including the junction's energy spectrum, effective capacitance and I — V characteristics. 

The corresponding critical behavior differs, in a number of important aspects, from the previously studied cases of 
a superconducting junction with the local Josephson energy as well as that of a normal juncion. However, in a general 
agreement with the earlier analyses of these systems the presence of a sufficiently strong quasiparticle and Cooper 
pair tunneling turns out to be instrumental for observing the Josephson behavior (albeit, possibly, not at the longest 
time scales). 

This research was supported by NSF under Grant DMR-0349881 and ARO under Contract DAAD19-02-1-0049. 
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FIG. 1: Phase diagram of the anisotropic XY-model. Renormalization group trajectories are shown by arrows, and the dashed 
(dotted) line represents the phase boundary between insulating and (super)conducting phases in the case of Ohmic (sub-Ohmic) 
dissipation. 



